Axial dipolar dynamo action in the Taylor-Green vortex 
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We present a numerical study of the magnetic field generated by the Taylor- Green vortex. We show 
that periodic boundary conditions can be used to mimic realistic boundary conditions by prescribing 
the symmetries of the velocity and magnetic fields. This gives insight in some problems of central 
interest for dynamos: the possible effect of velocity fluctuations on the dynamo threshold, the role 
of boundary conditions on the threshold and on the geometry of the magnetic field generated by 
dynamo action. In particular, we show that an axial dipolar dynamo similar to the one observed in 
a recent experiment can be obtained with an appropriate choice of the symmetries of the magnetic 
field. The nonlinear saturation is studied and a simple model explaining the magnetic Prandtl 
number dependence of the super/sub critical nature of the dynamo transition is given. 

PACS numbers: 47.65.-d, 47.20.Ky, 52.65.Kj 



I. INTRODUCTION 

The generation of magnetic field by the flow of an elec- 
trically conducting fluid has been mostly studied to un- 
derstand the magnetic fields of planets and stars [Ij. The 
increase in computing power has allowed the study of 
these fluid dynamos in almost realistic three-dimensional 
geometries although far from the parameter range rele- 
vant for astrophysical or geophysical flows or even labo- 
ratory experiments [2j. 

Three successful laboratory experiments have been 
performed so far: fluid dynamos have been observed with 
an helical Ponomarenko-type flow [3j or with a array of 
helical flows with the same helicity [4J. In these two first 
experiments, the flow lines were strongly constrained by 
the boundaries, thus restricting turbulent fluctuations to 
small scales, and the observed dynamos were found in 
very good agreement with the ones obtained by the action 
of the mean flow alone. A third example of fluid dynamo 
has been provided in the laboratory by the VKS exper- 
iment, i.e., a von Karman (VK) flow of liquid sodium 
[5J. In contrast to the previous dynamo experiments, 
the geometry of the generated magnetic field strongly 
differs from the predictions made by kinematic dynamo 
codes using the mean flow alone. The mean magnetic 
field has to leading order a dipolar structure with its axis 
along the axis of rotation of the propellers that drive 
the flow [6j, whereas kinematic dynamo codes using the 
mean flow alone predict an equatorial dipole. Thus, the 
VKS dynamo is not generated by the mean flow alone. 
It has been proposed that non axisymmetric fluctuations 
related to the propellers driving the flow play a crucial 
role in the generation of an axial dipolar magnetic field [7J 
and this mechanism has been illustrated by a kinematic 
dynamo simulation using a non axisymmetric model flow 
in a cylindrical geometry [8j. 

Two different choices can be made for the numerical 
simulations of dynamos. One possibility is to consider 



a flow confined in a finite domain and to use boundary 
conditions as realistic as possible. Another one is to use 
periodic boundary conditions on the velocity and mag- 
netic fields. In the later case, the Taylor-Green (TG) flow 
has been widely studied. The TG vortex is a standard 
turbulent flow often used in numerical computations [9j 
that is closely related to the experimentally studied VK 
swirling flow [TDHISJ. The relation between the VK flow 
and the TG vortex is a similarity in overall geometry: a 
shear layer between two counter-rotating eddies. The TG 
vortex, however, is periodic with impermeable free- slip 
boundaries (present as mirror symmetries) while the ex- 
perimental flow takes place between two counter-rotating 
coaxial impellers and is confined inside a cylindrical con- 
tainer. The TG vortex also obeys a number of additional 
rotational symmetries. 

Dynamo action in the TG flow was studied some time 
ago [13] and it was found that the most unstable growing 
magnetic mode is equatorial (perpendicular to the axis 
of rotation) and breaks the additional rotational symme- 
tries. Thus enforcing all the geometric symmetries of the 
TG flow was found to be non-favorable to dynamo action. 
Since that time all numerical studies of dynamo action 
in the TG vortex were performed using general periodic 
codes (not enforcing any of the symmetries). The gen- 
erated magnetic field was always found to be equatorial 
[14-16J. 

The new idea that we explore in the present paper is 
to enforce the mirror symmetries of the TG flow (that 
confine the flow) but not enforce the rotational symme- 
tries of the magnetic field (that forbid the appearance 
of magnetic dipoles). With mirror symmetries enforced, 
the corresponding boundary conditions on the magnetic 
field can be related to either electrically insulating [34J 
or perfectly conducting boundaries [17j. In the present 
case (without the rotational symmetries) there are 2^ = 8 
possibilities for the magnetic boundaries conditions. 

The paper is organized as follows: after recalling the 
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MHD equations and the TG forcing, we present the sym- 
metries of the TG flow in section [TTl We show that there 
exists only one choice of mirror symmetry for the veloc- 
ity field that is compatible with the equations of motion. 
It corresponds to free-slip boundary conditions on each 
boundary of the cube [0,7r]^. In contrast, in each direc- 
tions two independent choices exist for the magnetic field, 
either corresponding to an electrically insulating bound- 
ary, or to a perfectly conducting boundary. We show that 
we can always restrict the kinematic dynamo study in a 
symmetric subspace and that there are six different cases 
to consider. In section |III| we first shortly describe the 
numerical procedures and explain how the growth rates 
of the different dynamo modes are determined. Then, 
we consider how the imposed symmetries on the veloc- 
ity field affect the dynamo threshold. We show that in 
the absence of constraints, the fiow can spontaneously 
break the TG forcing symmetries, thus generating tur- 
bulent fiuctuations that inhibit dynamo action although 
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involving larger kinetic energy. In section IV , we consider 
the different symmetry constraints for the magnetic field 
and their respective dynamo efficiency. We show that the 
lower threshold is obtained when lateral boundaries are of 
different nature (two conducting and two insulating). In 
addition, the lowest threshold is obtained when the upper 
and lower ones are insulating. In this case, the threshold 
value and the neutral mode geometry correspond to the 
ones already reported in a general periodic code without 
symmetry constraints [13j. We observe that the geometry 
of the dynamo modes strongly depends on the symmetry 
constraints that are chosen. With all boundaries elec- 
trically insulating, we obtain the generation of an axial 
dipolar magnetic field, as the one observed in the VKS 
experiment. The nonlinear saturation of this axial dipole 
is studied in section |V| where we also present a simple 
model explaining the super/sub critical nature of the dy- 
namo transition when the magnetic Prandtl number is 
varied. Discussion and conclusions are given in section 

ED 



II. THEORETICAL PRELIMINARIES 

A. MHD equations and Taylor- Green forcing 

The magnetohydrodynamics (MHD) equations for an 
unit density fiuid and an incompressible fiow read in 
terms of the velocity v and magnetic induction b (in 
units of Alfven velocity): 

StV + v-Vv = -VP+j X b + i/Av + f (1) 
dth = V X (v X b) + T^Ab (2) 

together with V-v = = V-b;j = Vxbis the current 
density, P is the pressure, f is the mechanical forcing, u is 
the kinematic viscosity and 77 is the magnetic diffusivity. 



(3) 
(4) 



(where Ey{t) is the kinetic energy, Eij{t) is the magnetic 
energy and ( ) stands for spatial average over the do- 
main). The total energy is conserved by these equations 
in the ideal case {v = r] = and f = 0). 

Considering a flow that is 27r-periodic in all spatial 
dimensions, the kinematic Reynolds number Re and the 
magnetic Reynolds number Rcm are deflned as 



Re = 



LVr 



Rem = 



LVr 



V 



(5) 



where the root-mean square velocity is Vrms 
and the characteristic length L is defined by 



L = 2tt 



(6) 



where the kinetic energy spectrum Eyik^t) (such 
that Eyit) = ^]^Ey{k^t)) is obtained by summing 
^|u(k',t)p on the spherical shells A:- 1/2 < |k'| < A:+l/2 
(u(k) is the Fourier transform of the velocity). Anal- 
ogously, the magnetic energy spectrum is denoted by 
Ei){k^t) and verifles Ei)(t) = Ei){k^t). The magnetic 
Prandtl number of the fluid is given by P^ = u/rj = 
Rcm/Re. 

We now turn to the deflnition of the external driving 
volumic force f in ([T]) that balances the energy dissipation 
and allows to reach a statistically stationary state that 
is needed to sustain dynamo action. Following reference 
[13] , we force the system by setting in ([T]) 

i = f{tW, (7) 
where v^^ is the Taylor-Green [18j initial data 



TG 



(sin(x) cos{y) cos(z), — cos(x) sin(y) cos(2:), 0), 

(8) 

and f{t) is determined by imposing that the projection 
of V on the mode v^^ is flxed at all times to its initial 
value. 

In the Navier-Stokes problem (Eq. ([T]), with j = b = 0) 
is well-known [9j that a number of the symmetries of v^^ 
are compatible with the equations of motion. They are, 
flrst, rotational symmetries of angle tt around the axis 
{x = z = 7r/2) and {y = z = 7r/2); and of angle 7r/2 
around the axis {x = y = 7r/2). Another set of symme- 
tries corresponds to planes of mirror symmetry: x = 0, tt, 

= 0, TT and z = 0, TT that form the sides of the so-called 
impermeable box which conflnes the flow. However, it 
is important to note that the TG symmetries of the so- 
lution can be spontaneously broken^ in the sense that a 
small non-symmetric component of the initial data can 
grow and eventually break the symmetry of the solution 
|T3| . The symmetry breaking of the TG flow will be nu- 



merically studied in Sec III B 
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B. Mirror symmetries revisited 

The main new idea of the present paper is to enforce 
the mirror symmetries of the TG flow (that confine the 
fiow) but not enforce the rotational symmetries of the 
magnetic field (that forbid the appearance of magnetic 
dipoles). With mirror symmetries enforced, the corre- 
sponding boundary conditions on the magnetic field can 
be related to either electrically insulating or perfectly 
conducting boundaries |T7] (see the discussion at the end 
of Sec. II D). In the present case (without the rotational 
symmetries) it will turn out that there are new (mixed) 
possibilities for the magnetic boundaries conditions. We 
now proceed to explore all the TG mirror symmetries 
that are dynamically compatible with the full MHD equa- 
tions ( p!pl ). 

The MHD equations ([l][2| are invariant under refiec- 
tions with respect to a plane. We now turn to consid- 
erations pertaining to this propriety. Let us first de- 
fine the standard refiection transformation of the vector 
r = (ri,r2,r3) = {x^y^z) about the plane = by 
with 



S^(ri,r2,r3) = (-ri,r2,r3), 
S^(ri,r2,r3) = (ri,-r2,r3), 
S^(ri,r2,r3) = (ri,r2,-r3). 



(9) 
(10) 
(11) 



Note that is its own inverse. 

The action of the refiection operation on a vector field 
h(r) is defined by 

R"(h(r)) = S"h(S"r) (12) 

which explicitly reads, in the case of the the z = plane: 



^ha:{x,y,zy 

I hy{x,y,z) 
.hz{x,y,z) 



hx{x,y,-z) 
hy{x,y,-z) 
-hz{x,y,-z). 



(13) 



We will say that the vector field h has a mirror-symmetry 
with respect to the plane if it is even (R"h = — h) or 
odd (R"h = h) under transformation (12). Note that if 



a vector field is even (odd) with respect to a plane, then it 
is perpendicular (parallel) to that plane. When a vector 
field that is odd (or even) with respect to simultaneous 
transformations with respect to all planes, it is said that 
it transforms as a vector (or pseudo -vector). Note that 
the Taylor- Green vortex is odd with respect to all planes: 

j^a^TG ^ ^TG a. 

The following properties are true for any fields g and 



We now turn to study the mirror symmetries of the 
MHD equations. Applying to ([l][2| and writing v • 
Vv = (V X v) X V + V(^v^), it directly follows from 
Eqs. dMjlSliel that v = R"v and b = R"b satisfy 



d^-y + (V X v) X V 



dth 



1 



-V(P 
V X (v X b) 



^v2) + (Vb)xb 
R"f 



7?Ab, 



(17) 
(18) 



with P(r) = P(S"r). Therefore the only symmetry of the 
velocity that is compatible with the MHD equations ([l]- 
|2| is R"f = f and v = v. This can be easily understood 
with the following simple geometrical argument based on 
vortex rings that are known to have a self-induced motion 
along their axis. Suppose that there exist a symmetry 
plane containing a vortex ring, for such a configuration 
the vorticity is parallel to the plane. The related velocity 
field is thus an even field with respect to that plane. It 
is immediately noticed that the symmetry will not be 
preserved by the Navier-Stokes equation: the vortex ring 
will leave the plane under it's self-induced motion thus 
breaking the symmetry. Another possibility is to place 
the vortex ring perpendicular to the plane, making the 
velocity parallel, therefore odd with respect to that plane. 
The self-induced motion will then respect the symmetry. 

Note however that the magnetic field must be sym- 
metric but it can be either even or odd . Therefore, for 
a mirror symmetric velocity field, the magnetic field has 
two possible symmetries for each plane that are com- 
patible with the MHD equations. We will see below that 
the two possibilities correspond to insulating (I) and per- 
fectly conducting (C) magnetic boundary conditions (of 
the free-slip type) for the even and odd cases respec- 
tively. Note that, as the the fields are defined in a 27r- 
periodical box, the planes x = ir^y = tt and z = tt are 
also mirror symmetry planes with the same symmetries 
than X = O^y = and z = 0. In the following, depending 
on the symmetries of the magnetic field, we will refer to 
these planes as walls of type I (even) or C (odd). With 
this definition, the symmetry planes of the Taylor- Green 
vortex are of type CCC (vorticity is perpendicular to the 
box [0,7r]^). 



C. Projection on symmetric magnetic fields 

Let us now define the projectors into symmetric func- 
tions with respect to the planes = and = tt. 



R"(V X h) = -V X (R"h), 
R"(g X h) = -(R-g) X (R-h), 
R^Ag = A(R-g). 



(14) 
(15) 
(16) 



From Eq.(14) it is directly observed that if h is even 



with respect to a plane then V x h is odd with respect 
to that plane and vice versa. 



(19) 



Q? = -(/-sR"), 



where a stands for x^y^z^ R" is defined by Eq.([T2| and 
s = ±1. Note that by construction R"(Q^h) = -sQJh 
and therefore Q%.h. is an even vector (if 5 = +1) or odd 
vector {s = —1) with respect to the planes r^ = and 

Try = TT. 
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As the above discussion about the two possible mag- 
netic mirror symmetries was done for each plane Tq, = 
independently, we now define the projector over the most 
general mirror-symmetric magnetic fields corresponding 
to a different choice for each plane. 



(20) 



with s = (sx^Sy^Sz) G { — 1,1}^. Note that by construc- 
tion 



p(-l -1 -l)^TG ^ ^TG^ 



(21) 



Let us now consider a completely mirror-symmetric 27r- 
periodical velocity field. With this choice v is confined 
into the impermeable box [0, tt]^: there is no fiow crossing 
the boundaries. This mimics the rigid boundary condi- 
tions of experiments that are also impermeable. We now 
explain how the different kinds of magnetic symmetries 
are related to perfectly conducting and insulating bound- 
ary conditions. If the magnetic field is even with respect 
to one of the mirror-symmetry planes, we call this plane 
a insulating (I) wall because the current j = V x b is 
parallel to (or contained in) the wall [17J. Analogously 
if b is odd the wall is called conducting (the current is 
perpendicular to (or crosses) the wall. 

We now show that any 27r-periodical magnetic field can 
be decomposed in a finite sum of mixed symmetrical- 
periodical vector and pseudo-vector fields that corre- 
spond respectively to insulating and perfectly conducting 
magnetic boundary conditions. Thus the complete study 
of the infiuence of symmetries (boundary conditions) on 
the TG dynamo threshold is reduced to a finite number 
of possible cases. 

Let Q be the space of 27r-periodic functions, and = 
P^Q the projection of ft over the subspaces with mirror 
symmetries given by s. For example if b G then 
j = V X b is normal to the planes z = and z = tt and 
parallel to the other ones. Therefore is the 27r- 

periodical function space such that the planes x = 0, x = 
TT^y = O^y = TT are insulating and the planes z = 0, z = tt 
are conducting. 

It can be demonstrated by straightforward (but te- 
dious) computations that 



pSipS2 



I. 



(22) 
(23) 



We thus have that Q decomposes as the direct sum 

n = 0n,-, (24) 

S 

and therefore the general-periodic function space is 
then decomposed in eight mirror-symmetric periodical- 
function spaces. For the sake of clarity, in the following 
we will denote each space by its respective type of 
walls. For instance the case is labeled IIC. 



D. Decomposition of kinematic dynamo on 
symmetry classes 



The kinematic dynamo is the linear instability of the 
MHD equations linearized about a given velocity field Vg. 
For a stationary Vg, the growth rate of the magnetic 
field is given by the eigenvalue problem Ch = a^h with 



rb = V X (vg X b) + T^Ab. 



(25) 



Observe that the eigenvalue depends on r] and thus 
is a function of Rcm (see Eq.([5])). The dynamo thresh- 
old is defined by the critical magnetic Reynolds number 
Re^^* such that cr5(Re^^*) = 0. Note that this number is 
not necessary unique due to windows of instability [15]. 
In this work we refer to Re^^* as the smallest critical 
magnetic Reynolds number. 

As a consequence of the symmetry invariance of the 
MHD Eqs.([l][2|, it is straightforfard to show that the 
symmetry projectors Pg* commute with jC. Therefore if bo 
is an eigenvector associated to the eigenvalue <j^, by the 
decomposition of the identity (23), there exist one sym- 



metry So such that P^-^bo is also a eigenvector associated 
to <j^. We thus have for the threshold of general-periodic 
dynamo 



Re^^^* =min{Re^^^*'}, 



(26) 



where Re^^** is the critical magnetic Reynolds number of 
the linear problem (25) restricted to Qg. In other words. 



a bifurcating mode of the magnetic field has a certain 
mirror-symmetry. It is thus natural to study the dynamo 
bifurcation for each symmetry by restricting ( 25 ) to each 



Q^. This can be easily done if each mirror-symmetry is 
imposed to the magnetic field by applying the projectors 

Since the forcing ^ is invariant by rotation of 7r/2 
along the x = y = 7r/2 axis, the spaces ICI and ICC are 
respectively equivalent to CII and CIC. There are thus 
only six independent cases to study. The symmetries III, 
CCC, IIC, CCI, ICI and ICC wih be individually studied 
in SecHVl 



Note that we are dealing throughout this paper with 
periodic fields, with physical boundaries being replaced 
by symmetry planes. In contrast, in the case of a phys- 
ical fiuid confined into a finite domain a magnetic field 
at an insulating boundary should be matched with a po- 
tential field outside the fiow domain. The case that we 
call "insulating" also mimics a boundary with a magnetic 
permeability much larger than that of the fiuid, i.e. liq- 
uid sodium inside an iron vessel. In that case matching 
with a potential field is not required and the magnetic 
field is perpendicular to the boundaries. 
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III. 



A. 



SYMMETRIES OF THE VELOCITY FIELD 
AND DYNAMO THRESHOLD 

Numerical procedures and determination of 
growth rates 



a)io° 



Numerical solutions of Eqs.([l|2| are efficiently pro- 
duced using the pseudospectral general-periodic code 
GHOST [19j, that is dealiased by spherical spectral trun- 
cation using the 2/3-rule [20j. Thus a run with resolution 
has a maximum wavenumber /Cmax = Resolu- 
tions used in this works vary from 64^ to 256^. The equa- 
tions are evolved in time using a second-order Runge- 
Kutta method, and the code is fully parallelized with the 
message passing interface MPI library. We implemented 
into GHOST both the constant velocity forcing ^ and 
the projectors (20). 

The TG vortex ([8| is used as initial data for Eq.Q, 
eventually adding a small non-symmetric random part 
when studying symmetry breaking (see below section 
niB). A small, spectrally band-limited random seed of 
given magnetic energy is used as initial data for Eq.([2|. 

To correctly resolve the MHD Eqs.([l][2| spectrally, a 
fast decay at large k (faster than algebraic) of the en- 
ergy spectrum is required. This condition (called spec- 
tral convergence) is quantitatively determined by fitting 
the exponential decay of the energy spectra Ey{k^t) and 
E}){k^t) by a law of the form Ce~'^^^ that amounts to a 
simple Lin-Log linear regression. The value of Skma^ fur- 
nishes a measure of spectral convergence. For instance. 
Fig. [T] shows a numerical simulation where the magnetic 
field is well resolved, as apparent on the Lin-Log plot of 
E}){k). A fit of the data in Figjlja leads to S^kmax > 10. 
In all the runs presented in this work, we always have 
^^max > 2.5, with this condition we ensure that the fields 
are well resolved and that there is no spurious numerical 
effect on the observed dynamo instabilities. 

Once the correct resolution of MHD equations is en- 
sured the next step is to observe the behavior of the mag- 
netic energy Ei)(t). Within the linear theory of kinetic 
dynamo Ch = a^b (see Eq.([25|)) we expect an exponen- 
tial behavior of Ei){t) with a growth rate a^. For each 
run the magnetic energy is thus fitted with a law of the 
form E})(t) = Ce^^^ to determine a^. Figure [ijb shows a 
typical temporal evolution of Ei^it) and the correspond- 
ing fit gives the value = 0.017. The data in Fig jl] thus 
correspond to a typical dynamo action run. 



B. Symmetry breakings of velocity field 



As already stated at the end of Sec |II A[ the TG sym- 
metries of the solution to Eq.Q can be spontaneously 
broken, in the sense that a small non-symmetric compo- 
nent of the initial data can grow to order one and com- 
pletely break the symmetry of the solution. Two types 
of symmetry breaking will be considered in this section: 
ffist confinement breaking and second (with confinement 
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FIG. 1: a) Magnetic energy spectrum of a dynamo run: case 
III with Re = 30, Rcm = 80. b) Corresponding temporal evo- 
lution of magnetic energy (|3|. The fits used to determine 
Gh = 0.017 and 6kmax = 10.72 (for the present simulation) 
are displayed as (dashed) straight lines of Figs, a) and b). 



enforced) the breaking of the rotational symmetry of an- 
gle 7r/2 around the axis {x = y = 7r/2). We now turn to 
a numerical study of these points. 



i. Spontaneous confinement breaking 

To study the spontaneous breaking of the mirror sym- 
metries of the TG flow (that confine the flow) we compare 
runs performed with and without projecting the velocity 
field after each timestep using ( [21] ). The time dependence 
of the kinetic energy at Re = 30 is shown in|2]a, where it 
is seen that (when non enforced) the mirror symmetry is 
broken with an increase of the kinetic energy Ej^{t) due 
to fluctuations of the velocity field. Mirror (confinement) 
symmetry breaking is seen to take place around Re = 10 
in Fig|2jb, where the dependence of the kinetic energy 
(time-averaged over statistically stationary values) on Re 
is displayed. 



2. Spontaneous 7v/2 rotation breaking 

We now turn to the numerical study of another bifur- 
cation of the steady-state velocity field that takes place 
when confinement ( [21] ) is enforced. Stable and unsta- 
ble steady states are followed by making use of Newton's 
method, the linear equations being solved by the stabi- 
lized bi-congugate gradient algorithm [21j, see [22j and 
references therein. The result of such computations is 
displayed in Fig[3ja, where it is apparent that a pitch- 
fork bifurcation is present at critical Reynolds number 
Re^ = 22. Physical 3D visualizations of the velocity 
fields are made using; VAPOR [35J and presented in Fig[3] 
b-c. By computing the difference between the stable and 
unstable branches obtained near the bifurcation by New- 
ton's method the bifurcating mode can be numerically 
obtained. The dominant Fourier component v^^ (PF 
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FIG. 2: (Color online) a) Temporal evolution of kinetic en- 
ergy with (red) and without (blue) confinement (21) imposed 
at Re = 30. b) Dependence of kinetic energy on the Reynolds 
number for the symmetric solid (red) line and non-symmetric 
dashed (blue) line runs (time-averaged over statistically sta- 
tionary values). 
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FIG. 4: (Color online) Dynamo growth rates a correspond- 
ing to the symmetric and non symmetric velocity fields of 
Fig|2^ as function of the magnetic Reynolds number Rem- 
Blue crosses: v and b non symmetric, Red stars: v symmet- 
ric and b non symmetric. 



standing for pitchfork) is found in this way to be 



PF 



sin(x) cos(?/) cos(z) 
cos(x) sin(y) cos(z) 
-2 cos(x) cos(?/) sin(z) , 



(27) 



It is straightforward to check that v^^ is odd under 7r/2 
rotation around the axis {x = y = /2). As the equations 
and the non-bifurcated states are both invariant by this 
transformation the bifurcation is a pitchfork (see Figjsja) 
that breaks the 7r/2 rotational invariance around the axis 
{x = y = 7r/2). 
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FIG. 3: (Color online) a) Bifurcation diagram: kinetic en- 
ergy H/v (k = 2) as function of Re. A pitchfork bifurcation 
is clearly present {Ev{k = 2) is quadratic in the bifurcat- 
ing mode v^^). b-c) Visualization of Taylor- green station- 
ary states at Re = 30: non-bifurcated (b) and bifurcated 
(c) velocity fields. Streamlines are colored accordingly to the 
magnitude of the velocity field, varying form magenta (low 
intensity zones) to yellow (high intensity zones). 



C. Effect of velocity fluctuations on magnetic 
energy growth rates 

The symmetries of the confined velocity fields reduce 
the fluctuations of the velocity field, as observed in Figj2j 
Using the methods outlined in section III A to determine 
the magnetic growth rates, we observe that this reduction 
enhances the magnetic instability and reduces the criti- 
cal magnetic Reynolds number. This is apparent in Figj4] 
where the growth rates for the non-symmetric and sym- 
metric velocity field are compared as functions of Rcm at 
fixed kinematic Reynolds number Re = 30. 



IV. EFFECT OF MAGNETIC BOUNDARY 
CONDITIONS 



A. Magnetic boundary conditions and dynamo 
threshold 



We now focus on a confined velocity field (with mir- 
ror symmetries imposed) and study the influence of the 
different kinds of magnetic boundary conditions on the 
dynamo threshold. A number of runs have been per- 
formed at kinematic Reynolds numbers Re = 30, Re = 
150 and magnetic Reynolds numbers varying between 
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Re^ 

^ FIG. 6: (Color online) Dynamo growth rates cr^ (symmetric 

velocity field) obtained with a non-symmetric magnetic field 
FIG. 5: (Color online) Dynamo growth rates at (symmetric compared with the CIC magnetic symmetric case. Kinetic 
velocity field) as function of the magnetic Reynolds number Reynolds number Re = 30 (a) and Re = 150 (b). 
Rem for the 6 possible symmetries of the magnetic field at 
kinetic Reynolds number Re — 30 (a) and Re = 150 (b). 



Rem = 10 — 300 with resolutions of 64^ — 256^ grid points. 
For each run, the growth rate is measured and the whole 
ensemble of data is presented in Fig. |5j This is our main 
quantitative result on growth rates. 

By performing a linear interpolation of we obtain 
the critical magnetic Reynolds number Re^^* that are 
given in Table [T] for Re = 30. Observe that critical 



Case 


ICI 


ICC 


lie 


III 


CCC 


CCI 




9 


26 


66 


73 


231 


254 



TABLE I: Critical magnetic Reynolds number for the different 
walls and symmetric and no n- symmetric cases at Re = 30. 
Values obtained by linear fit of at. 

magnetic Reynolds numbers vary from Re^^* 10 to 
ReS^* - 250 for the different kind of walls. 

The most favorable cases correspond to mixed 
insulating-conducting lateral walls (ICI and ICC). These 
configurations allow for a magnetic field crossing the box 
in the direction perpendicular to the insulating walls and 
current crossing in the other direction. The less favor- 
ables cases turn out to correspond to the lateral perfectly 
conducting walls (CCC and CCI). Note that, see the dis- 
cussion at the end of Sec |IID[ the condition we call insu- 
lating corresponds to a magnetic field perpendicular to 
the boundary. This case can be achieved experimentally 
by using a ferromagnetic boundary. 

Fig. [5] also displays the growth rate for Re = 150. At 
this Reynolds number, the velocity field is turbulent and 
the fluctuations increase the dynamo threshold. Observe 
that the case ICI is still the most unstable. 



Note that, as a consequence of the direct sum decompo- 
sition presented in Sec II C| (Eq.(24))), the most unstable 
case (ICI) and the case with no symmetries imposed on 
the magnetic field should have the same critical magnetic 
Reynolds number. This condition is indeed verified (see 
Fig. [g]), both cases having Re^^* = 9. Furthermore, the 
property also implies that the growth rates of these 
two cases are equal. This is apparent in Figj6]a, where 
the two growth rates are plotted for Re = 30. Both curves 
are almost identical for all magnetic Reynolds numbers 
used in this work. The slight difference at Rcm = 50 
may be due to imprecisions in the numerical determi- 
nation of small growth rates caused by interferences be- 
tween competing modes. At Re = 150 the curves only 
coincide qualitatively. At this relatively high kinematic 
Reynolds number the flow is turbulent and the instan- 
taneous growth rate (defined as in Fig. [l]b) fluctuates 
thus explaining the slight discrepancy between the two 
curves. 



Magnetic boundary conditions and geometry of 
unstable modes 



We now turn to study the geometry of the unstable 
modes near the dynamo threshold by generating 3D vi- 
sualizations of the magnetic field and current. In all the 
figures presented in this section the velocity field is sym- 
metric. A resolution of 256^ is used and visualization are 
made using VAPOR. In all visualizations we represent 
(color online) the magnetic field lines (in red), the cur- 
rent (in yellow) and a density plot of the highest magnetic 
energy zones. The axis of forcing (z-axis) is indicated by 
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a blue arrow and the equatorial directions {x and y axes) 
by red and green arrows located on one corner of the box. 

The two cases corresponding to a non-symmetric mag- 
netic field (a) and ICI walls (b) are first compared in 
Figj?! Observe that, for both cases, the magnetic field is 




FIG. 7: (Color online) 3D visualizations (magnetic field in 
red, current in yellow and density plot of highest magnetic 
energy zones) of the growing modes: (a) Case v symmetric 
and b non-symmetric, Rcm = 10. (b) ICI case, Rcm = 10. 



very similar: the magnetic field lines are mainly in one 
equatorial direction. The only difference being their re- 
spective orientation, with respect to the x and y axes. 
For the ICI case the magnetic field must cross the in- 
sulating walls, imposing its direction, while for the non- 
symmetric case this direction is chosen randomly by the 
fiow between the cases ICI and CIC (recall that the cases 
ICI and CIC are equivalent by 7r/2 rotation). This field 
is very similar to the one reported in [13j. 

The visualization of cases ICC, IIC, CCC and CCI are 
displayed in FigjSj Note that both ICC and IIC cases 
contain axial magnetic fields but the magnetic field lines 
do not cross the top and the bottom of the box, a pos- 
sibility that is not allowed by the boundary conditions 
(C walls on top and bottom). The cases CCC and CCI 
present complicated structures. 

Fig. [9)a displays another of the main results of this 
paper: an axial dipole that is obtained with III walls. 
The magnetic lines are clearly oriented along the 2:-axis 
(see the zoom view displayed on[9jb). To the best of our 
knowledge this is the first time that an axial dipole has 
been observed using periodical boundary conditions. 




FIG. 8: (Color online) 3D visualizations (magnetic field in 
red, current in yellow and density plot of highest magnetic 
energy zones) of the growing modes: (a) ICC case, Rcm = 30 
(b) IIC case, Re^ = 80, (c) CCC case, Re^ = 300, (d) CCI 
case. Rem — 300. 




FIG. 9: (Color online) 3D visualizations (magnetic field in 
red, current in yellow and density plot of highest magnetic 
energy zones) of the growing modes: (a) III case, Rcm = 80. 
(b) Zoom of figure (a) at the center of the box. 



V. NON-LINEAR SATURATION OF AXIAL 
DIPOLAR MAGNETIC FIELD 

We have studied so far the kinematic dynamo problem 
and have found how the first unstable magnetic mode, 
when Rcm is increased, depends on the boundary con- 
ditions. We present in this section the different station- 
ary regimes that are obtained when the growing unstable 
mode saturates due to the back reaction of the Lorentz 
force on the velocity field. This study is restricted to 
the case III for which an axial dipole is the first unstable 



magnetic mode. 

The critical Re^^* bifurcation line, observed in nu- 
merical dynamo s imulatio ns (III case) around the stable 
branch of section IIIB2 (see Fig. [sja), is presented in 
Fig. 10 a (lower curve (6 and 2) with diamonds). The 
vertical line at Re^ = 22 (1 and 5) corresponds to the 
hydrodynamic pitchfork bifurcation discussed in section 
nil B 21 These two bifurcation lines intersect at a so-called 
codimension-2 bifurcation point. Dynamical regimes in 
the vicinity of this point are displayed as the phase por- 



traits of Fig 10 b that will be discussed below. 




FIG. 10: (Color online) a) Qualitative drawing of the bifurca- 
tion lines. In the vicinity of the codimension-two point where 
these lines intersect, they are given by model (29). Diamonds 
correspond to data coming from numerical simulations (blue 
lines), lines without markers are qualitatively drawn (green 
lines). Each bifurcation line is labeled with a number (see 
text) and each region of parameter space is labeled by I- VI . 
Red triangles corresponds to some simulations used to identify 
the different regions of the parameter space. Filled (empty) 
markers denote a non-vanisnhing (vanishing) magnetic field in 
saturated regimes. Half- filled triangles denote points in the 
bi-st able zone, b) Phase diagrams (in A-B space, see section 
VD) corresponding to regions I- VI of (a). 



A. Dynamo shutdown by velocity bifurcation 

Let us first emphasize that the nonlinearly saturated 
regime is not always related to the linearly growing mode. 
This type of behavior, obtained for Re = 30, is displayed 
in Fig|TT]a. For Re^i = 80, i.e. above the bifurcation 
threshold Re^,'^* 



73 (region III in Fig 10 a), the mag- 
netic energy displays a growing phase (after some tran- 
sient) but then decreases and vanishes in the long time 
limit. Figure [TT]b shows that this change of behavior of 
the magnetic energy occurs concomitantly with a varia- 
tion of the kinetic energy of the flow. 

We have checked that this variation is related to the 
transition from the symmetric flow taken as initial condi- 
tions to the one that breaks the 7r/2 rotational invariance 
around the axis {x = y = 7r/2) (see section IIIB2). Al- 
though the symmetric flow is hydrodynamically unstable 
for Re = 30, the initial conditions are such that it has 
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FIG. 11: Case III: Temporal evolution of kinetic and magnetic 
energy at Re = 30 and Rcm = 80. Vanishing of saturated 
magnetic field is clearly observed at large times. 



not yet broken its symmetry when the magnetic field be- 
gins to grow. It is the Lorentz force that drives the flow 
to one of its bifurcated states with broken 7r/2 rotational 
invariance. The bifurcated flow having no dynamo capa- 
bility for Rcm = 80, the magnetic field then decays. Note 
that a similar phenomenon was observed for the magnetic 
field generated by a flow in a spherical domain |23| . 



B. Supercritical dynamo 

The complex transient behavior of the magnetic field 
being related to the interaction with the hydrodynamic 
pitchfork bifurcation reported in section [IIIB 2[ we now 
study the dynamo bifurcation at Re = 20, below this hy- 

a). 



drodynamic bifurcation threshold (region VI in Fig 10 



As shown in Fig. 12 a the linear unstable magnetic mode 
saturates at finite amplitude thus displaying a supercriti- 



cal pitchfork bifurcation. The kinetic energy on Fig. 12 b 
is slightly reduced by the saturation mechanism. It is also 
apparent on Fig. [T2]c-d that the magnetic field lines are 
slightly more twisted in the saturated regime, as can be 
checked visually by focusing on the two magnetic field 
lines that are colored in green in the figure. 

We now consider how the magnetic energy Q depends 
on the fluid parameters above the dynamo threshold. 
Dimensional arguments imply (b^) = (v^) /(Re, Rcm) 
where / is an unknown function. Close to a supercrit- 
ical bifurcation threshold, we expect that (b^) depends 
linearly on Rcm — Ke^^*, thus 



(b2)^(ve2)5(Re) 



r^e r^e 

Re"'' 



(28) 



where (vc^) is related to the kinetic energy density at 
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FIG. 12: (Color online) Case III: Temporal evolution of mag- 
netic (a) and kinetic energy (b) at Re = 20 and Rem = 80. 
c) 3D visualization of magnetic field at linear growth phase 
t = 400. d) Visualization of saturated magnetic field at 
t = 3000. Two of the magnetic field lines have been colored 
in green in order to emphasize the twisting of the saturated 
field. 



bifurcation threshold. It is expected that ^(Re) tends 
to a constant in the limit of large Re and is inversely 
proportional to Re for small Re [24] . Using v^^^ = 2Ey / 3 
to estimate (vc^), we find that ^(Re) ~ 3 increases by 
about 20% when Re is varied from 18 to 21 (data not 
shown). The high Re number scaling thus is observed 
even at moderate values of Re (compared to experiments 
with liquid metals). 



C. Bistability 

For Re > 22, a bistable region can be found in which, 
depending on the initial conditions, we can get both a 
dynamo and purely hydrodynamic regimes. For instance, 
this case was observed (data not shown) at Re = 26 and 
Rem = 100 (region IV in FigjlOja). 

Starting from a dynamo in such a bistable regime, we 
followed a line at Rcm = 100 by increasing the Reynolds 
number up to Re = 120 without loosing the dynamo 
(data not shown). 

A hysteresis was also observed, for Re > 22, by varying 
Rcm at fixed Re (data not shown). 



D. Codimension-2 bifurcation model 

We next present a simple explanation for the super/sub 
critical nature of the dynamo transition as the kinetic 
Reynolds number (or equivalently the magnetic Prandtl 
number of the fluid) is varied. In our simulations, this 
phenomenon is strongly related to the presence of a pitch- 
fork bifurcation of the flow for Re = Re^ = 22. 

For Re = Re^ and Re^ = Re^^*(Rec) = Re^, we have 
a codimension-two bifurcation. In its vicinity, the dy- 
namo and the hydrodynamic instabilities compete, thus 
generating various dynamical regimes. 

Denoting by A{t) the real amplitude of the bifurcating 
velocity field v^^ and B{t) the real amplitude of the bi- 
furcating magnetic field b, we write coupled amplitude 
equations for A and B in the vicinity of the codimension- 
two bifurcation. The form of these equations is con- 
strained by symmetry requirements, A —A (pitchfork 
bifurcation of the velocity field) and B —B (b — b 
symmetry of the MHD equations). 

Keeping the nonlinear terms to leading order, we get 



A = XA-aAB^ -A^, 
B = jiB - pA^B - B^. 



(29) 



The coefficients of the cubic nonlinearities have been 
taken negative in order to get supercritical pitchfork bi- 
furcations for the hydrodynamic instability in the ab- 
sence of magnetic field {B = 0) and for the dynamo in- 
stability when Re < Re^ and thus A = 0. The modulus 
of these coefficient can be taken equal to 1 by appropriate 
scalings of the amplitudes A and B. A and fi are functions 
of Re and Rem that vanishes at the codimension-two bi- 
furcation point {Re^^Re^). To leading order, we have 
X (X Re — Re"^ and ja oc Rem — Re^'^^. 

The fixed points of the system (29) are 
(0,0), (±a/A, 0), (0, ±y7^) and the mixed modes 
{±y^{X - a/i)/(l - a/3), ±^{ii - pX)/(l - a/3). The 
different types of bifurcation diagrams have been studied 
in detail p5] . 

The cases of interest are presented in Fig{To]b. They 
correspond to a/3 > 1 (in order to have unstable mixed 
modes) with a and /3 positive (in order to prevent the ex- 
istence of subcritical mixed modes for A and /i negative) . 
Then the globally stable solution (0, 0) for A < and 
/i < undergoes a supercritical pitchfork bifurcation to 
(±>/A, 0) (respectively (0, ±y7^)) for A = (respectively 
/i = 0). 

The corresponding bifurcation lines are labelled (1) 
(respectively (6) in Fig. 10). For A > (respectively 



/i > 0), the second pure mode bifurcates for /i = (re- 
spectively A = 0) but is unstable (bifurcation lines (2) 
and (5)). 

The mixed modes exist between the lines A = a/i and 
(1 = /3X (corresponding to (3) and (4) in Fig. 10). The 
key point is that they are unstable with respect to the 
pure modes. The system is thus bistable in this parame- 
ter range. 
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The bifurcated hydro dynamic regime as weh as the 
dynamo state are both hnearly stable. When Rcm is 
increased for Re > Re^, the hydrodynamic regime bi- 
furcates to a dynamo state on hue (4). If Rcm is then 
decreased from this state, the dynamo is suppressed on 
hue (3), thus displaying an hysteresis. 

An hysteresis can be also observed by varying Re for 
Rcm > Re^. If the fluid velocity is increased at con- 
stant magnetic Prandtl number, i.e. if one follow a line 
Rcm = PmRe in parameter space, we expect a super- 
critical dynamo bifurcation for Pm large enough and a 
subcritical one for Pm small. 



Let us finally remark that the model ( 29 ) is expected to 



be valid only in a neighborhood of the codimension-two 
point (i^e^, Re^). Indeed many other secondary bifurca- 
tions can take place away from this point. Because of this 
and also due to the limited number of runs that were per- 



formed, the lines in Fig 10 a are only qualitatively drawn. 



VI. DISCUSSION AND CONCLUSION 

It has been often claimed that too many symmetries 
of the velocity and/or magnetic fields inhibit dynamo 
action. This claim probably results from several anti- 
dynamo theorems that have been found since Cowling 
who showed that an axisymmetric velocity field cannot 
generate an axisymmetric magnetic field (for a review, 
see [26j). Thus, the magnetic field should break axisym- 
metry when generated through dynamo action by an ax- 
isymmetric velocity field. Another class of anti-dynamo 
theorems is even more restrictive and forbids the dynamo 
action of some velocity fields (for instance planar flows, 
i.e. velocity fields with only two non zero cartesian com- 
ponents) whatever the geometry of the generated mag- 
netic field. 

We have shown here that symmetries do not always in- 
hibit dynamo action, but, on the contrary can sometimes 
enhance it: 

- symmetry constraints on the velocity field can lead 
to lower dynamo threshold by inhibiting the development 
of hydrodynamic instabilities and related turbulent fluc- 
tuations that sometimes reduce the efficiency of dynamo 
action, 

- symmetry constraints on the magnetic field can lead 
to an unchanged dynamo threshold provided they are 
chosen in the appropriate manner. 

It has been shown that the dynamo threshold Re^^* of 
the TG flow, increases when the kinetic Reynolds num- 
ber of the flow is increased on some intermediate range, 
and this has been related to the development of turbulent 
fluctuations U4l[T5l[27]. Using symmetry constraints on 
the velocity field shows that phenomenon without ambi- 
guity: the level of turbulent fluctuations is lower with 
the symmetric velocity field and, correspondingly, the 
dynamo threshold is lower too, although the mean ki- 
netic energy is larger for the velocity field without sym- 
metries. This shows in a simple and clear-cut way that 



velocity fluctuations inhibit dynamo action by the TG 
flow. The same phenomenon has been shown analyti- 
cally for the dynamo generated by a fluctuating Roberts 
flow [28]. However, it should be kept in mind that a time- 
dependent velocity can also generate new dynamo modes 
that do not exist in the absence of fluctuations. 

We have also shown that symmetry constraints on the 
magnetic field can be used to mimic realistic bound- 
ary conditions in the framework of numerical simula- 
tions with periodic codes. It was found that the dy- 
namo threshold and the geometry of the growing mag- 
netic mode strongly depend on the choice of symmetry 
constraints, and thus on the related boundary conditions. 
The lowest dynamo threshold was obtained with lateral 
boundaries of different nature, allowing a magnetic field 
to cross the box perpendicularly to the current (ICI case) . 
This case, where the dominant component is an equato- 
rial dipole, strongly reminds the geometry of the mag- 
netic field that was numerically generated by using the 
mean flow measured in a VK geometry with counter- 
rotating propellers [29", "30]. This emphasizes the sim- 
ilarity between VK and TG flows. Note however that 
the mean flow component in the VK geometry is axisym- 
metric and cannot drive an axisymmetric magnetic field 
because of Cowling's theorem. A magnetic field with an 
equatorial dipolar component provides a simple way to 
break axisymmetry. In contrast, the TG forcing is not 
axisymmetric and thus does not enforce such a strong 
constraint on the generated magnetic field. Indeed, when 
symmetry constraints related to infinite magnetic per- 
meability are implemented for all boundaries (III case), 
the generated magnetic field involves a dominant com- 
ponent which is an axial dipole. This corresponds to 
the geometry of the magnetic field generated in the VKS 
experiment, where it was ascribed to the presence of non- 
axisymmetric velocity fluctuations (T). 

Direct numerical simulations of dynamos generated by 
an axisymmetric s2 — t2 forcing in a sphere, a config- 
uration similar to the VK forcing, have shown that an 
equatorial dipole is observed at low Re, for which the 
flow is axisymmetric, whereas an axial dipole is obtained 
at higher Re, i.e. in the presence of non-axisymmetric ve- 
locity fluctuations [31j. We observe that the TG forcing 
provides a different scenario: the axial dipole does not 
seem to be favored by turbulent fluctuations (compare 
Fig. [5] a and b). Thus, slight deviations from axisym- 
metry, present in the TG forcing, are enough to generate 
a magnetic field with a dominant axial dipolar compo- 
nent, even at small Re, provided appropriate boundary 
conditions are simulated. 

Finally, we have studied the nonlinear saturation of 
the axial dipolar magnetic field (III case). We have ob- 
served that the dynamo bifurcation can be subcritical in 
some parameter range. We have shown that this phe- 
nomenon can be explained by the presence of an hydro- 
dynamic instability that competes with the growing mag- 
netic mode. This is not unlikely for Pm of order 1 which 
is the case of most numerical simulations and provides a 
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general mechanism for subcritical dynamos. It would be 
of interest to check whether the subcritical dynamo bi- 
furcations reported previously in non confined TG flows 
[32j, flows forced in a sphere [23] or simulations of the 
geodynamo [33j, are also related to a similar mechanism. 
Other open questions concern the dynamics of the mag- 
netic field above the dynamo threshold. Breaking the 



symmetry of the TG flow in order to generate reversals 
of the axial dipole or other dynamics resulting from non- 
linear coupling between different magnetic modes deserve 
to be studied. 
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